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O \ I have studied the variabihty properties and locahzation of 3EG J1746-2851 based 

Q \ on EGRET data of the observing periods 1-4. Using corrections for know systematic 

' problems and performing various consistency checks I find no evidence of variabihty 

■ with an amphtude exceeding 30% with the possible exception of viewing period 429, 
for which a strong soft excess is observed. 3EG J1746-2851 is displaced from the exact 

^ . Galactic Center towards positive Galactic longitudes. Sgr A*, the center of Sgr A 

2 . East, the pulsar J1747-2958, and the TeV 7-ray source observed with HESS seem to be 

\Q \ excluded as possible counterparts at the > 95% level. 
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1. Introduction 

The Galactic Center region contains a many potential sources of high-energy 7-ray emission 
(Melia &: Falcke 2001). A point source of GeV-scale 7-ray emission, 3EG J1746-2851, has been 



\ detected with EGRET (Mayer-Hasselwander et al. 1998), but no firm identification with sources 

^ ' known in other frequency bands has been made to date. Many models have been published to ex- 

plain 3EG J1746-2851, e.g. invoking an advection-dominated accretion flow onto the supermassive 
^ ' black hole known as Sgr A* (Mahadevan et al. 1997; Markoff, Melia and Sarcevic 1997), the decay 

■ of neutralinos (Berezinsky et al. 1992), the acceleration of quasi-monoenergetic electrons in the 

radio arc (Pohl 1997), or a pulsar located on the line-of-sight to the Galactic Center (McLaughlin 
& Cordes 2003). 

Very high-energy 7-rays have also been detected from the Galactic Center region with the 
Whipple 10m telescope (Kosack et al. 2004), the CANGAROO telescope (Tsuchiya et al. 2004), 
and the HESS array (Aharonian et al. 2004). The HESS data indicate a power-law spectrum 
with photon index a = 2.21 it 0.09, an extrapolation of which does not fit the GeV-scale flux 
measured with EGRET. The source of TeV-scale 7-ray emission shows no evidence of variability 
and is coincident with Sgr A* to within about 0.02°. The lack of a cut-off in the 7-ray spectrum 
seems to disfavor models involving the annihilation of dark-matter particles (Bergstrom et al. 1998; 
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Ellis et al. 2002; Gnedin &: Primack 2004), for a very high mass > 12 TeV of the hypothetical 
neutralino would be required. 

More conventional models for the high-energy emission from the Galactic Center region have 
also been proposed. Electron acceleration by plasma turbulence near the sTipcrmassive black hole 
(Schodcl ct al. 2002) may be a very efficient process (Liu et al. 2003). Atoyan & Dermer (2004) 
argue that the combination of ADAF-type accretion with a subrelativistic MHD outflow can explain 
the TeV-scale emission as well as that at X-rays and lower frequencies, but 3EG J1746-2851 would 
be a separate source. Fattuzzo & Melia (2003) discuss Sgr A East, a nonthermal radio source 
with a supernova-like morphology located near the Galactic center, as a possible counterpart of 
3EG J1746-2851, and Crocker et al. (2004) attempt to relate 3EG J1746-2851 , the TeV-scale 7-ray 
source, and anisotropies in the UHE cosmic-ray intensity to Sgr A East. 

One of the fundamental questions is whether or not the sources seen in the GeV and TeV bands 
are identical. The analysis of EGRET data for the Galactic Center region is known to be affected by 
systematic problems, some of which have been understood since the detection of 3EG J1746-2851 
was announced in 1998. In this paper I re-analyze the available EGRET data with a view to infer 
the variability properties and the localization as accurately as reasonably possible. For that purpose 
I introduce corrections in the way the diffuse foreground model is treated in the likelihood analysis, 
which are outlined in Sec. 2. I do not change the foreground model per se. Various consistency 
checks are performed using data for different energy bands and from different observing dates. 

2. Modeling the galactic foreground emission 

I have first analyzed the combined data sets of the EGRET observing periods 1-4 in an area 
of ±60° in latitude and longitude around the Galactic Center. Using all 80 sources, that are listed 
in the Third EGRET Catalogue (Hartman et al. 1999) to be located in an area of ±40° in latitude 
and longitude around the Galactic Center, I have run the LIKE software (Mattox et al. 1996) to 
simultaneously determine the multipliers Guuit Gbiss for the standard model of galactic diffuse 
7-ray emission (Hunter et al. 1997) and isotropic emission, respectively, and the fluxes of these 80 
sources while keeping their positions fixed. The multipliers Guuit and Geias were allowed to vary 
freely. The purpose of the exercise was two-fold: 

• to test whether or not the combination of the eighty 3EG sources and the foreground model 
is an acceptable fit to the entire data set. If there are 7-ray excesses with test statistic 
TS > 16 within 30° of the Galactic Center, we may need to incorporate them in the analysis 
of individual viewing periods. 

• to derive an appropriate range of parameter values for the foreground model. The galactic 
diffuse emission cannot be variable, and therefore the same model, i.e. the same values for 
GmuH and Geias that we obtain from analyzing the combined data set, should be used in 
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the likelihood analysis of individual viewing periods. However, Guuit a-^d Geias must be 
allowed to vary within the range of statistical uncertainty of their best-fit values and within 
the uncertainty range of the spark-chamber efficiency correction (Esposito et al. 1999). 

In contrast to the procedure used to derive the Third EGRET Catalogue, I did not include below- 
threshold excesses with TS > 9 in the analysis. The source fluxes obtained for the combined data 
set of the EGRET observing periods 1-4 are generally very similar to those listed in Hartman et al. 
(1999). The most prominent exception is 3EG J1744-3011, located at I = -1.15° and b = -0.52°, 
with 

3 EG S = 9Aa F{> 100 MeV) = (63.9 ± 7.1) 10"^ ph./cmVsec 
This analysis S = 7.6a F{> 100 MeV) = (50.5 ± 6.9) 10"^ ph./cmVsec (1) 

For the Galactic Center source, 3EG J1746-2851, the flux and source significance determined with 
the combined data set are virtually identical to those listed in the catalogue. 

3 EG 5=17.5 0- F(> 100 MeV) = (119.9 ±7.4) 10~%h./cmVsec 

This analysis S = 17.5 a F{> 100 MeV) = (118.6 ± 7.3) 10"^ ph./cmVsec (2) 

At the position of the Galactic Center source the best-fit parameters of the foreground model are 

GMuit = 0.948 GBias = 3.405 (3) 

In the entire area of ±20° in latitude and longitude around the Galactic center these parameters 
have values not exceeding the range [0.93,1.18] for GMult and [2., 4.] for Geias- The high value 
for the multiplier of isotropic emission, GBiasj is possibly related to inverse-Compton emission not 
accounted for in the standard model of galactic diffuse 7-ray emission. 

The highest intensity of diffuse galactic 7-rays is observed in direction of the Galactic Center. 
Any systematic errors associated with a mismodeling of the galactic emission will therefore be 
most severe in this region of sky. Systematic effects will also be more important at low 7-ray 
energies than at high energies on account of the strong energy dependence of EGRET's point-spread 
function (PSF). Given the average photon flux above 100 McV (Eq.2) and the similarity between 
the apparent 7-ray spectrum of 3EG J1746-2851 (Merck et al. 1996) and that of the diffuse emission, 
we can infer that the mean measured intensity (or count density) associated with 3EG J1746-2851 
within the solid angle element containing 68% of the photons, J7(68%), is approximately 

/gc(68 %) ~ 0.68 ^^^(g8%)''^^ - ^"^ ■ Ph-/cmVsec/sr (4) 

which is small compared with the mean intensity of diffuse emission in that solid angle element. In 
the energy band above 1 GeV the mean intensity associated with 3EG J1746-2851 is actually more 
than a factor of two higher than that of the diffuse emission, solely on account of the narrower PSF. 
We must therefore espect that the high-energy band above 1 GeV is relatively clean of systematic 
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problems, as far as the determination of flux and position of 3EG J1746-2851 are concerned, and 
that the low-energy bands and the very wide bands, e.g. > 100 MeV, are subject to possibly serious 
systematic uncertainties. 

One such systematic uncertainty arises from a simplification in the EGRET standard data 
analysis chain. There one uses allsky maps of modeled galactic diffuse 7-ray emission, in which 
the effect of the PSF is already incorporated. These maps are then multiplied with the exposure 
map (the time integral of the effective area) of the viewing period in question to produce a map of 
the expected event distribution that can be compared with the data. So the PSF is applied before 
the exposure, whereas it should be the other way around. Unpublished earlier Monte-Carlo studies 
by the author have indicated that this simplification may be problematic for 3EG J1746-2851, if 
the pointing direction is in the Galactic Plane at galactic longitude / ~ ±20°, so that the Galactic 
Plane is radially oriented in the EGRET field-of-view, and the effective area of the instrument is 
rapidly decreasing with off-axis angle at the location of the point source in question. 

For each viewing period I have created customized maps of diffuse galactic 7-ray emission, 
that are corrected for this simplification. Using a simple analytical representation of the spectrum 
of the diffuse galactic 7-ray emission. 



as well as the effective area distribution, the energy dispersion function, and the tabulated PSF 
data from the pre-launch calibration files for the dominant EGRET observing mode 74, I have 
calculated a mean PSF for each energy interval, for which data analysis would be performed. The 
mean PSF thus derived was then applied to the expected true count distribution, i.e. the product 
of the model of galactic emission, GALDIF, and the exposure map, EX, to derive the expected 
count distribution in the data. 



where a and (3 denote the true angular coordinates of the incoming 7-rays, whereas 9 and (f) are the 
measured angular coordinates of the 7-rays. The standard maps of diffuse galactic 7-ray emission, 
that can be obtained from the CGRO Science Support Center, are calculated assuming a single 
power-law spectrum with index 2.1, thus implying a PSF of different width that used in my analysis. 



Table 1 displays a list of all viewing periods considered in this study. For each viewing pe- 
riod the data have been analyzed in the energy bands 100-300 MeV, 300-1.000 MeV, >100 MeV, 
>300 MeV, and > 1.000 MeV. I have also constructed a combined data set for all viewing periods 
with a view to derive a detailed spectrum and to pinpoint the localization of the source. 




for E <l GeV 
for E>1 GeV 



(5) 
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Table 1. A list of all viewing periods considered in this study. The table also shows the pointing 
direction in galactic coordinates and the observing time. 
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3.1. Consistency checks and the spectrum 

The spectra of EGRET sources are usually determined on the basis of a likelihood analysis in 
ten standard energy ranges. Table 2 lists the integrated photon flux and the v flux in the ten 
narrow energy bands and five wide energy bands based on a likelihood analysis of the combined 
data set. The spectrum depends very weakly on the assumed position of the source, if that is 
varied by bl < 0.1° in the Galactic Plane. Figure 1 shows the v Fy spectrum above 70 MeV. The 
parameters of the diffuse foreground model, GMult ^-nd Ceias) were allowed to vary freely. Table 2 
also lists the integrated photon flux for all energy bands that one obtains using the standard maps 
of diffuse foreground emission. The combined data set has an exposure distribution that peaks 
within 2° of the Galactic Center and is fairly symmetric. The standard maps and the modified 
maps of galactic diffuse emission differ only marginally. 

The photon flux for the energy bands below 100 MeV is strongly dependent on the particulars of 
the likelihood run, e.g. on the analysis parameter RANAL, and thus discredited. An anticorrelation 
between the multiplier GmuIi of the diffuse emission model and the photon flux attributed to 
3EG J1746-2851 can be observed, which indicates a similarity of the expected angular distributions 
of events for the diffuse emission and 3EG J1746-2851, respectively. The results of the likelihood 
analysis at energies below 100 MeV do therefore not merit our consideration in the variability and 
localization study. 

One also notes that the integrated flux above 100 MeV is significantly higher than derived in 
the allsky analysis (Eq.2). The only difference between the two likelihood analyses was the inclusion 
in the allsky analysis of all viewing periods not listed in table 1, which leads to a significantly more 
homogeneous distribution of exposure across the sky. However, the exposure distribution for the 
combined data set as listed in table 1 is wider than that of any individual viewing period. So if 
the likelihood analysis of the combined data set is apparently inaccurate at low 7-ray energies for 
which the PSF is wide, there is no reason to assume it would fare better for any individual viewing 
period. 

Another consistency check would be a comparison of the integrated flux obtained for a wide 
energy band with the sum of the fluxes in the corresponding narrow energy bands. While the 
high-energy bands are consistent with each other, we note that for the wide band 100-10.000 MeV 
the flux is about twice the sum of fluxes obtained in the corresponding narrow energy bands. The 
flux determined for the 300-10.000 MeV band is about 14% higher than the sum of the fluxes for 
the corresponding narrow energy bands. The chance probability of that deviation is about 15%, so 
given the number of trials we have to expect one deviation of that amplitude on statistical grounds. 
Significant problems can thus only be noted for the energy band 100-10.000 MeV. 



- 7- 



Table 2. The spectrum of 3EG J1746-2851 as derived using the combined data set. The 
integrated photon flux in each energy band is given in units of 10~^ cm~^s^^. The u F^, flux is in 
units of IQ-^ MeV cm ^ s ^ and derived assuming a photon spectrum cc E ^ . The results of the 
hkchhood analysis at energies below 70 McV should be considered as upper limits for systematic 
reasons. Upper limits are given for a 2a confidence level. For comparison the fourth column gives 
the integrated photon flux derived using the standard maps of diffuse emission. 
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Fig. 1. — The v Fi, spectrum of 3EG J1746-2851 as derived from an analysis of the combined data 
as hsted in table 1. The spectrum is similar to that of the galactic diffuse emission. 
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3.2. Search for variability 



In the preceding subsection we have reported on a number of inconsistencies in the results of 
a Hkelihood analysis of the combined data set, which indicate systematic problems. The affected 
data in the energy bands 100-10.000 MeV and below 100 MeV will therefore be excluded from the 
following analysis. 

We now want to analyze the data for individual viewing periods with a view to infer the 
variability properties of 3EG J1746-2851. Earlier publications (e.g. Nolan et al. 2003) reported 

some evidence for variability based on a likelihood analysis of EGRET data in the energy band 
100-10.000 MeV that we found affected by systematic problems. Here we use only data for energy 
bands, for which we see no clear evidence for systematic errors. 

In the search for variability we therefore analyze EGRET data in the energy bands 100-300 
MeV, 300-1.000 MeV, 1.000-10.000 MeV, and 300-10.000 MeV. In the narrow energy bands that 
we have used to derive the spectrum of 3EG J1746-2851 the statistical significance of detection is so 
low that an analysis of data of individual viewing periods does not provide meaningful results. We 
keep the position of the source fixed, for that cannot vary. Our analysis of the combined data set 
yielded / = 0.11° and b = -0.04° as best-fit position of 3EG J1746-2851, which we henceforth use 
as fixed location in the search for variability. Small changes in the assumed position by 61 < 0.1° 
have virtually no impact on the results of the variability analysis. 

The intensity of diffuse emission should be constant over time, and so should be the values 
of its scaling factors in the likelihood analysis, Guuit and GBias- The uncertainty of the spark- 
chamber efficiency correction introduces some variations in the best-fit values of the scaling factors, 
that should have an amplitude around 10% (Esposito et al. 1999). For each data set we ran the 
likelihood analysis twice, first with the scaling factors fixed to the values found in the allsky analysis 
(Eq. 3), and then allowing the scaling factor to freely vary. Variable scaling factors should provide 
a better representation of the diffuse galactic emission, for they can balance errors in the spark- 
chamber efficiency correction and they can correct inappropriate large-scale structure in the model 
of diffuse emission. On the other hand, one may expect an anti-correlation between the best-fit 
flux from 3EG J1746-2851 and the best-fit value of Guuit, which is the dominant parameter for the 
intensity of the expected diffuse 7-ray emission from the inner galaxy. Therefore the mean relative 
amplitude of variations in GmuU over N viewing periods 



should not significantly exceed 10%. We thus monitor the variations in the best-fit values of 
CMuit to ensure that variations in the measured flux from 3EG J 1746-2851 are not caused by an 
unrealistically extreme value of GmuH- 




N 



(7) 
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3.2.1. The energy band 300-1 0. 000 Me V 

The results of the hkelihood analysis runs for the individual viewing periods as well as for the 
combined data set are summarized in table 3. The scaling parameter of diffuse emission, GmuH) 
has best-fit values that are characterized by 

G^t = 1.021 aGunit = 0.106 (8) 

so the scatter in GyiuXt is at a level commensurate with the systematic error in the absolute deter- 
mination of egret's effective area at the time of measurement, thus not indicating any additional 
systematic problem. 

We can perform a x^-tcst to determine how well the measured fluxes in the 24 viewing periods 
are compatible with a constant flux. The best-fit constant flux would be 

F{> 300 MeV) = 45.74 . 10"® cm'^ s"^ (9) 

for which the minimal value of the x^-sum and the chance probability of drawing the measured 
distribution given a constant flux, the goodness-of-fit, would be 

= 30.0 Pehance = 15% (10) 

This indicates some scatter in the distribution of flux values beyond what must be expected on sta- 
tistical grounds. We can estimate the amplitude of that scatter by Gaussian adding of a systematic 
uncertainty 5 = y F %o the flux measurement error. The amplitude of scatter is then approximately 
given by the value of y, for which the minimal value of the x^-sum equals the number of degrees of 
freedom (22 in this test). The resulting estimate for the amplitude of variations is 

ybest fit = 0.197 (11) 
which is just twice the expected level of absolute calibration error. 



3.2.2. The energy hand 300-1.000 MeV 

This energy band is part of the energy band, the results for which we have just discussed. It 
is thus not statistically independent. The scaling parameter of diffuse emission, GmuMi has best-flt 
values that are characterized by 



GMuit = 0.875 crGMuit = 0.118 (12) 

so the scatter in Gyiwit is again at a level commensurate with the systematic error in the absolute 
calibration. 

The best-flt constant flux would be 

F(300 - 1.000 MeV) = 25.08 . 10"® cm^^ g-^ (13) 



- 11 - 



Tabic 3. Results of the likelihood analysis in the energy band > 300 McV. The position of 
3EG J1746-2851 was held fixed at I = 0.11 deg and b = -0.04 deg. The integrated photon flux F 
is given in units of 10~® cm~^s~^. Upper limits correspond to 95% confidence or 2a. 





Gm 


uit fixed 1 




GmvlH variable 


VP 


a 


F 


a 


F 




r 




K o 

o.z 


OQ n_i_^ o 


7 Q 

i .6 


44.9±7.0 


0.889 


Id 


o.y 


oO.Oihii. 1 


7 Q 
1 .6 


76.4±12.3 


1.129 




^ 
z.o 


^Sl /I4-07 A 


1 Q 

i.y 


49.0±28.6 


1.129 




i.D 


oo Q-i-on 1 
zy.oitzu. 1 


i.O 


28.8±21.3 


1.059 


01 Q 


n 7 

U. 1 


^1 OK 7 


U.O 


<136.6 


0.988 


OOQ 


Q 

z.o 


/1 1 04_on Q 

4i.o±zu.y 


1 Q 

i.y 


37.7±22.0 


1.117 


OOi^ 
ZZD 


i.y 


1 n n_i_ 11 ^ 
iy.yitii.o 


o o 

z.z 


25.9±12.6 


0.969 




O 1 

Z.L 


QO Q_l_ 1 V 


O 1 

z.i 


37.1±19.4 


1.015 


OQ1 


O /I 


^^_i_on n 


/I 
Z.4 


62.6±31.2 


0.980 


OQO 
ZoZ 


O.D 


f;i Q-Ui Q n 


0.4 


64.5±14.0 


1.029 


3023 


2.4 


29.7±13.4 


1.9 


24.7±14.3 


1.068 


323 


4.9 


45.9±10.7 


3.9 


39.4±11.3 


1.122 


324 


3.6 


54.0±17.5 


3.2 


52.0±18.5 


1.069 


330 


5.0 


118.2±29.6 


4.3 


110.5±31.1 


1.086 


332 


6.8 


69.1±12.1 


6.6 


72.2±12.9 


1.022 


334 


2.0 


27.7±14.9 


1.9 


27.5±16.0 


1.053 


3365 


2.4 


55.2±26.9 


1.9 


46.3±27.9 


1.120 


421 


1.8 


20.5±12.7 


2.5 


32.2±14.5 


0.938 


422 


1.2 


12.9±11.2 


3.2 


38.1±13.5 


0.822 


423 


1.9 


17.3±9.8 


3.7 


38.1±11.6 


0.897 


4235 


0.3 


< 37.2 


2.2 


34.2±18.0 


0.773 


429 


3.2 


56.7±20.7 


4.0 


77.8±23.3 


0.883 


508 


5.6 


88.7±19.7 


4.0 


69.5±20.2 


1.235 


625 


4.2 


43.1±11.7 


3.2 


36.4±12.6 


1.107 


all 


17.9 


44.7zb2.9 


17.8 


48.5±3.1 


1.011 



- 12 - 



for which the minimal value of the x^-sum and the chance probability of drawing the measured 
distribution given a constant flux, the goodness-of-fit, would be 

= 31.0 Pchance = 12-3% , (14) 

very similar to the results obtained for the wide energy band 300-10.000 MeV. 
The resulting estimate for the amplitude of variations is 

2/best fit = 0.346 , (15) 
which is nearly twice as much as in the wide energy band. 



3.2.3. The energy hand 1.000-10.000 MeV 

At these very high energies we expect the model of diffuse galactic emission to play a small 
role in the likelihood analysis. On the other hand, the small number of photons will reduce the 
statistical accuracy of the results. 

The scaling parameter of diffuse emission, Giviuit, has best-fit values that are characterized by 

G^t = 1.46 (tCmuH = 0.154 (16) 

so the scatter in GMuit is slightly larger than the systematic error in the absolute calibration. The 
large value of GmuH reflects the GeV excess in the diffuse galactic 7-ray emission (Hunter et al. 
1997). 

The best-fit constant flux would be 

F(> 1.000 MeV) = 12.9 • 10"^ cm'^ s"^ (17) 

for which the minimal value of the x^-sum and the chance probability of drawing the measured 
distribution given a constant flux, the goodness-of-fit, would be 

X^in = 25.3 Pehance = 33.5% , (18) 

indicating no evidence for variability at energies higher than 1 GeV. 
The resulting estimate for the amplitude of variations is 

ybest fit = 0.182 , (19) 

which here is more an estimate for the sensitivity of the method used in the search for variability, 
for additional fiuctuations witha an amplitude much higher than this value would be required to 
significantly increase the minimum x^-sum in Eq. 18. 
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3.2.4- The energy range 100-300 MeV 

This data in this energy range should contain a large number of 7-ray events, but are potentially 
problematic because of the large width of the PSF. Systematic problems were evident at lower 
energies, and we do not know very well to what extent the data at 100-300 MeV are affected. 

The scaling parameter of diffuse emission, Gnuit; has best-fit values that are characterized by 

G^t = 0.897 aGuuit = 0.132 (20) 

so the scatter in GmuW is marginally larger than the systematic error in the absolute calibration. 

The best-fit constant flux would be 

F(100 - 300 MeV) = 29.9 • 10~^ cm^^ s~^ (21) 

for which the minimal value of the sum and the chance probability of drawing the measured 
distribution given a constant flux, the goodness-of-fit, would be 

xLn = 43.58 Pchance = 0.6% . (22) 

The resulting estimate for the amplitude of variations is 

ybest fit = 0.71 . (23) 

A close inspection of the flux distribution reveals that all this evidence of variability is caused by 
one data point, that for viewing period 429. If we neglected that measurement, our results would 
be 

F(excl. VP429) = 26.1 • 10"^ cm^^ s'^ (24) 

and 

xLn = 21.92 Pchance = 46.5% cxcl. VP429, (25) 

i.e. no evidence for any variability. So the question is whether or not the measurement of a high 
flux during VP429 is realistic. 

The measured flux from 3EG J1746-2851 during VP429 is 

i^VP429(100 - 300 MeV) = (242.6 ± 46.1) • 10"^ cm^^ s"^ (26) 

What was the flux in the other energy bands? In all cases it was around 60% above average with 
the statistical significance of the excess being around la. Equation 26 indicates that in the energy 
band 100-300 MeV the flux was a factor of 8 above average with a statistical significance of 4.6(7. 
If the outburst were real, the source must have an soft flare spectrum, as shown in figure 2. 

What was the best-fit value of the scaling factor of diffuse emission, GMuit, for VP429? It is 
G'Muit(VP429) = 0.66, actually the lowest value of all viewing periods, which mandates a view at 
the results of the likelihood analysis for fixed Gjviuitj 

i^VP429(GMuit = 0.948) = (104.1 ± 39.9) • 10"^ cm"^ s"^ (27) 

which is still more than three times the average flux with a statistical significance of 2a. 
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Fig. 2. — The u Fi, spectrum of 3EG J1746-2851 during viewing period 429. Data shown with 
dotted or dashed error bars refer to wide energy bands and are not statistically independent of 
the data points plotted with solid error bars. The spectrum is much softer than the time-averaged 
spectrum shown in Fig. 1. 
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3.3. The localization of 3EG J1746-2851 

The likelihood ratio test, that is used to determine the flux and detection significance of point 
sources, can also provide an estimate of the source position and its uncertainty. For that purpose 
one calculates the point source test statistic 

TS = -2 [In Lo - In Li (a, /?)] (28) 

for a finely sampled array of test positions (cc, /?) (Mattox et al. 1996). Here Lq is the likelihood 

function for the null hypothesis and Li is the likelihood function assuming a point source at the 
position (a, (3). The increase in TS at the location of its maximum over TS at the true position is 
expected to follow a xi distribution, for the two angular coordinates are two additional degrees of 
freedom (Wilks 1938). Decrements from max(TS) by 2.3, 6.0, and 9.1 delineate the boundaries of 
regions in which the source is located with a confidence level of 68%, 95%, and 99%. 

Obviously the localization analysis can provide meaningful results only if max(TS) is not too 
small. We have therefore first analyzed the combined data set, which should provide the best 
statistical accuracy. As in our search for variability (Sec. 3.2), we exclude data at energies below 
100 MeV and in the wide energy band 100-10.000 MeV, because they appear to be affected by 
systematic problems. 

In the energy band 100-300 MeV 3EG J1746-2851 is detected with a significance of only 3.5a, 
too low to provide constraining results, so here we only report that the localization estimate for 
that energy band is compatible with those obtained using data at higher photon energies. 

In the figures 3 and 4 we show source localization confidence contours for the energy bands 
300-1.000 MeV and > 1.000 MeV, respectively, in comparison with the position of the central part 
of the galactic center arc and that of Sgr A*. To be noted from the figures is that in the two 
independent energy bands the likelihood analysis indicates a source position that is significantly 
displaced from the exact Galactic Center, or Sgr A*, or the position of the TeV 7-ray source 
(Aharonian et al. 2004), that coincides with Sgr A* to within 1.2 minutes of arc. 

The best-fit positions determined for the two energy bands are moderately well compatible 
with each other. However, both arc within the 68% confidence contour as determined for the 
100-300 MeV band. The 68% confidence contours for the energy bands 300-1.000 MeV and > 
1.000 MeV just touch each other, but there is a substantial overlap of the 95% confidence regions. 
The localization analysis of these two statistically independent data sets thus provides a much 
more consistent result than did the original study by the EGRET team (Mayer-Hasselwander et 
al. 1998), where even the 99% confidence contours for the 1-30 GeV and the 100-300 MeV bands 
didn't come close to each other. 

We can further test for systematic errors in two ways, which we discuss in turn. 
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Fig. 3. — Source localization confidence contours for the energy band 300-1.000 MeV. The square 
indicates the position of the central part of the Galactic Center arc and the asterisk denotes Sgr 
A*. The triangle and the diamond indicate the best-fit position in the energy bands 300-1.000 
MeV and > 1.000 MeV, respectively. 
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Fig. 4. — The results of the source localization analysis for the energy band > 1.000 MeV. The 
symbols indicate the positions of the same systems as in Fig. 3. 
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3.3.1. Consistency of localization in different energy bands 



The test statistic, TS, of two statistically independent data sets is additive. We can therefore 
co-add the TS fine maps for the energy bands 300-1.000 MeV and > 1.000 MeV to obtain a TS 
fine map for the total energy range > 300 MeV. The localization information thus derived should 
be compatible with the results of a likelihood analysis of the data in the band > 300 MeV, which is 
systematically different for two reasons. First it lacks the information which photons have a high 
energy and thus a narrow PSF, and therefore there is less information available. Second, the source 
spectrum may be different from the simple power law that is assumed to construct a band-averaged 
PSF, and therefore the PSF used in the likelihood analysis is the more likely to be inaccurate the 
wider the energy band. In Fig. 5 we compare the confidence contours derived from a likelihood 
analysis of the data in the > 300 MeV band with those obtained by adding the TS fine maps for the 
two narrower bands. There is a substantial overlap of the confidence regions, and the differences 
can be attributed to the displacement of the source localization at energies > 1.000 MeV relative 
to that at lower energies, an information that is preserved in the co-added TS fine map but not 
available to the likelihood analysis for the energy band > 300 MeV. We do not note any systematic 
problem that might arise from a possibly inaccurate PSF in the analysis of the wide energy band. 
Again, the position of Sgr A* and the position of the TeV 7-ray source seem to be excluded, and 
the GeV 7-ray source 3EG J1746-2851 appears to be located a fraction of a degree away from the 
exact Galactic Center. 



The localization accuracy obtained from the analysis of data for individual viewing periods will 
be worse than that derived using the combined data set, and the best-fit positions of 3EG J1746- 
2851 will be scattered over an area much larger than the confidence regions shown in Figs. 3-5. 
We can test whether or not the observed scatter in the best-fit source position is compatible with 
a common true position and the measured localization uncertainties. The confidence contours can 
usually be accurately fitted with an ellipse that is specified by a semimajor axis p, semiminor axis 
q, and position angle of the major axis cp (Thompson et al. 1995). The ellipse is centered on the 
centroid of the confidence location region, which thus serves as a position estimate. Mattox et al. 
(1997) noted that for strong sources the dependence of the logarithm of the likelihood of EGRET 
data upon the assumed source position closely follows a paraboloid, indicating a Gaussian error 
distribution and thus justifying the use of the x^-t6st. In contrast to Thompson et al. (1995) and 
Mattox et al. (1997) we use a fit to the 68% confidence contour. The displacement of the assumed 
(or measured) position from the true position is described by the distance r and the position angle 
9. Then the x^-sum is defined by 



3.3.2. Consistency of localization in individual viewing periods 




(29) 
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Fig. 5. — The results of the source localization analysis for the energy band > 300 MeV. The solid 
lines denote the confidence contours derived from a likelihood analysis of the data in the > 300 MeV 
band, whereas the dotted lines are obtained by adding the TS fine maps that result from likelihood 
analysis of the statistically independent data in the 300-1.000 MeV and > 1.000 MeV bands. The 
meaning of the symbols is the same as in Fig. 3. 
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A certain minimum in detection significance is necessary to provide a meaningful localization. Here 
we only consider data from viewing periods for which 3EG J1746-2851 was detected with at least 
VTS = 3. 

Table 4 lists the results of the location analysis in the energy range > 300 MeV for the 12 view- 
ing periods, for which the detection significance exceeded the required threshold. By minimizing 
the x^-sum as given in Eq. 29 we obtain for that energy range 

Xmin = 8.876 Ard.o.f. = 10 Pchance = 54.4% (30) 

The coordinates, for which the x^-sum is minimized, and their marginalized Icr-errors are 

I = 0.38° ± 0.09° b = -0.007° ± 0.075° (31) 

which is right inside the 68% localization contour derived from teh combined data set for the same 
energy band (Fig. 5). The scatter in the position estimates for the 12 viewing periods and the 
best-fit position thus determined are therefore entirely consistent with the statistical uncertainties. 
There is no evidence for an additional systematic error in the localization. 

4. Summary and discussion 

In this paper I have rc-analyzcd EGRET data of the observing periods 1-4 to derive the 
variability properties and localization of 3EG J 1746-2851. For that purpose I have introduced 
corrections in the way the diffuse foreground model is used treated in the likelihood analysis, and 
I have also performed various consistency checks using data for different energy bands and from 
different observing dates. Using only data for which I found no evidence of systematic problems, 
e.g. excluding the wide energy band 100-10.000 MeV and other low-energy bands, I obtained the 
following results: 

• There is little evidence of variability. The chance probability of drawing the measured dis- 
tribution of photon fluxes given a constant flux is always higher than 10%, corresponding 
to a x^/d.o.f. smaller than 1.35. The observed scatter in the measured flux values may 
correspond to variability at a level of 20-30%, if it were real. Thus I find an upper limit of 
< 30% to the relative amplitude of stochastic variability in the lightcurvc above 300 McV 
of 3EG J1746-2851. This conclusion would be corroborated by the data in the energy range 
100-300 MeV, if I neglected viewing period 429, for which one notes an apparently intense 
flare with a soft spectrum. I find some indications of systematic problems in the data of 
VP429 for the 100-300 MeV band, but no substantial indication that the observed high flux 
may be an artifact. 

• The data indicate that 3EG J1746-2851 is displaced from the exact Galactic Center towards 
positive Galactic longitudes, in qualitative agreement with the findings of Hooper & Dingus 
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Table 4. Results of the localization analysis in individual viewing periods and the energy band 

> 300 MeV. The position angle (j) = and (j) = tt/2 correspond to the direction of the 
galactic-longitude axis and the negative galactic-latitude axis, respectively. Wc only list viewing 
periods for which 3EG J1746-2851 was detected with at least VtS = 3. The best-fit positions are 
in degrees, the ellipse axes are in arcmin, and the major-axis position angle is in degrees. 
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(2002). This displacement is consistently found in maximum likelihood analyses of all energy 
bands above 300 MeV and also using data of individual viewing periods and the energy range 
300-10.000 MeV, in contrast to the contradictory results of Mayer-Hasselwander et al. (1998). 
The position of Sgr A*, the center of Sgr A East, the pulsar J1747-2958, and the position 
of the TeV 7-ray source seem to be incompatible at the > 95% level with the results of the 
localization analysis of 3EG J1746-2851. The Galactic Center arc or any other source located 
in the Galactic plane at longitudes 0.1° ^ / ^ 0.4° arc possible counterparts to the GeV 7-ray 
source. On account of the similarity of the spectra of 3EG J1746-2851 and that of the diffuse 
Galactic 7-ray emission it is possible, though very unlikely, that a compact, dense complex 
of interstellar gas, that is unaccounted for in the standard foreground model, is responsible 
for 3EG J1746-2851. The total mass of the gas cloud would have to be Mgas c± 5 • 10^ M© if 
located at the Galactic Center, and the cloud should have an apparent extent of Ogas ^ 0.3°, 
lest it doesn't appear point-like at high 7-ray energies. Such a massive gas complex should 
have been seen in the atomic fine-structure lines (Ojha et al. 2001; Rodriguez-Fernandez et 
al. 2004) or in CO lines (Kimc et al. 2002; Martin ct al. 2004). A somewhat smaller mass of 
unaccounted gas is required to shift the apparent position of a true point source, but even at 
the 20% level, i.e. Mgas — 10^ M0, the radio and infrared data do not indicate the existence 
of a previously unknown gas of the appropriate localization and compactness. 
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